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ABSTRACT 

Here we investigate the contribution of surface Alfven wave damping to the heating 
of the solar wind in minima conditions. These waves are present in regions of strong 
inhomogeneities in density or magnetic field (e. g., the border between open and closed 
magnetic field lines). Using a 3-dimensional Magnetohydrodynamics (MHD) model, we 
calculate the surface Alfven wave damping contribution between 1-4 Rq (solar radii), 
the region of interest for both acceleration and coronal heating. We consider waves with 
frequencies lower than those that are damped in the chromosphere and on the order of 
those dominating the heliosphere: 3 x 10~ 6 — 10 _1 Hz. In the region between open and 
closed field lines, within a few Rq of the surface, no other major source of damping has 
been suggested for the low frequency waves we consider here. This work is the first to 
study surface Alfven waves in a 3D environment without assuming a priori a geometry 
of field lines or magnetic and density profiles. We demonstrate that projection effects 
from the plane of the sky to 3D are significant in the calculation of field line expansion. 
We determine that waves with frequencies >2.8 xlCP 4 Hz are damped between 1-4 Rq. 
In quiet sun regions, surface Alfven waves are damped at further distances compared 
to active regions, thus carrying additional wave energy into the corona. We compare 
the surface Alfven wave contribution to the heating by a variable polytropic index and 
find that it an order of magnitude larger than needed for quiet sun regions. For active 
regions the contribution to the heating is tweny percent. As it has been argued that 
a variable gamma acts as turbulence, our results indicate that surface Alfven wave 
damping is comparable to turbulence in the lower corona. This damping mechanism 
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should be included self consistently as an energy driver for the wind in global MHD 
models. 

Subject headings: Sun: corona, Sun: magnetic fields, Sun: solar wind, waves 



1. Introduction 



The physical mechanisms behind the heating of the solar corona and the acceleration of the 
fast solar wind are two major unresolved issues in solar physics. Thermal heating alone is not 
sufficient to bring models in to agreement with observations of the lower corona and at Earth 
(jUsmanov fc Goldstein! 120031 ) . The acceleration of the sola r wind occurs pre dominantly within a 
few solar radii of the surface (jHartmann MacGregorlll980l ; lGrail et al.lll996l ). Additionally, Solar 
and Heliosqheric Observatory (SOH O ) observations hay e shown that ions are heated below 4 R Q 



(Kohl et al 



1998 



Esser et al.1 Il999l ). I Grail et al.1 (|l996l ) suggested that because the locations for 



the heating of the corona and the acceleration of the solar wind are the same, it is possible that 
the same mechanism could contribute to both. 

One source of heating is magnetic reconnectio n associated with flares and nanq nares at the solar 
surface. Using extreme ultraviolet observations, iPatsourakos h Klimchuld (120091 ) have suggested 
the heating of active regions is impulsive, and therefore could be associated with nanoflares. In the 
high speed solar wi nd, reconnection events current sheets and filaments have been suggested as a 
heating mechanism (jMattheaus et al.ll2003l ). 



In the lower corona, the region of interest in this work, damping of Alfven waves is known 
to produce nonthermal acceleration and to bring models into agreement with observations both 
near the Sun and at large distances ()Parkerlll965l ). The possible damping mechanisms fo r Alfven 



wave s in the photosphere to lower corona are n umerous and include nonlinear damp ing (jWentzel 



1989 



i ng (Hevvaerts Priestl 



1983; 



(Dc 



Parker 

Pontieu. Martens fc HudsonlhoOl 



2001 



Ofman k, Davilalll997r). t urbulent casc ade (|Hollweglll986l : iMattheaus et al.lll999r ). phase mix- 



1991) , Landau da mping (jHollwegl Il97ll ). neutral collisional 



Leake ct a 



), and surface Alfven wave damping (jlonson 



ion-cyclotron damping (llsenberg. Lee Hollwegj 



19781 ) (and references herein). 



Global magnetohydro dynamics (MHD) models that do not include waves employ different 
methods such as empirical heating functions or varied polytropic index distribution to drive and heat 
the solar wind and are able to match well with Ulysses, Yohkoh, Helios an d Advanced Composition 



Explorer data (jMikic et al.1 Il999l : iGroth et al 



2000 



Roussev et al.ll2003l ). To include Alfven wave 



damping in a realistic model of the solar environment, physical damping mechanisms must be 
specified. Some MHD model s use Alfve n waves without specifying a damping mec hanism. Global 
(jUsmanov &; Goldstein! 12000 ) and local (ICranmer. van Ballegooiien k, Edgarll2007l ) models include 
additional equations for the waves and prescribe an empirical damping length. These studies are 
benchmarked with Ulysses and Helios observations. Alternatively, one can include the wave energy 
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and momentum without damping (jLionello. Linker &; Mikid [2009j), and match SOHO Extreme 
Ultraviolet and Yohkoh soft X-ray observations. 

Of the numerous possible mechanisms for damping Alfven waves, those which we expect to be 
important for low frequency waves are nonlinear turbulent damping, phase mixing, surface Alfven 
waves. In this paper, we study the contribution of surface Alfven wave damp ing by utilizing a 3D 
global MHD simulation of a thermally-driven solar wind ([Cohen et al.l 120071 ). In this model, the 
lower corona is the inner boundary (chromosphere and transition region are not resolved.) Waves 
are not explicitly included in the model, nor can low frequency waves be resolved due to spatial 
and temporal resolution limits (to be discussed further in Section 2.2). 

We consider waves with frequencies lower than those that are completely damped in the chro- 
mosphere and on the order of those dominating the heliosphere: 3 x 10 -6 — 10 _1 Hz (periods 3 
seconds to 3 days). Alfven waves hav e been detected in the lower layers of the s olar atmosphere 



using both ground based observations (I Jess . et alj|2009f ) and the Hinode spacecraft (jOkamoto et al 
20071 : ICirtain et al.l [2QO7I : be Pontieu et al.ll2007h with periods 2-4 minutes (4.2 - 8.3 x 10" 3 Hz) 



Ground based observati ons indicate the presence of Alfven waves in the co rona with periods of five 
minutes (3.3 x 10" 3 Hz) (jTomczvk et alJboOTl ; rTomczyk & Mclntoshlbood ) . At 1 AU, the dominant 
wave power is m waves with periods of 1-3 hours (9.2 x 10 -5 to 2.8 x 10" 4 Hz) (Belcher fc Davisl 



1971 



)■ 



In this study we are interested in a damping mechanism that acts in the lower corona (| Grail et al 



19961 ). In the chromosphere, Alfven waves with frequencies above .6 Hz are damped by ion-neutral 



collis i onal damping and f requencies below 10 2 Hz were unaff ected ( De Pontieu. Martens Hudson 



200ll ; lLeake et al.ll2005l ). ICranmer &: van Ballegooiienl ( 2005) found that waves below 1 2 Hz were 



Cranmer &: van B allego oiien (120051 ) also 
Verdini fc Vellil (|200^ found 



not damped by any mechanism in the chromosphere, 
found that nonlinear damping occurred over the extended corona, 
waves 10 -6 — 10 -4 Hz were reflected by a gradient in the background Alfven profile, and dissipated 
not in the lower layers of the atmosphere, but over distance of a few solar radii. In the corona, 
the ion-cyclotron frequency is 10 4-6 Hz, so cyclotron resonance damping is not relevant for low 
frequen cy waves. Phase mixing of out going and reflected incoming Alfven waves has also been 
studied lsuz~uki Inutsukal (|2005l . bood ). 



Using a combination of t hree damping mechanisms (non linear damping, surface Alfven wave 
damping, and phase mixing), iJatenco-Pereira &: Opherl (|1989l ) were able to match observations of 
mass loss rates and terminal velocities for cool, giant stars. They applied their model to the Sun 
and were able to obtain coronal h eating and match wind velocity and Alfven wave power density 
observations in a ID simulation (jJatenco-Pereira. Opher Yamamotd 11994). In this paper, we 
extend their work on surface Alfven wave damping in a 3D simulation of the solar corona. 

Surface Alfven waves form on a magnetic interface - a finite thickness boundary separating 
two regions of plasma with a strong inhomogeneity in magnetic field and/or density. The Alfven 
wave in each region can interact, damp and transfer energy into the resonant layer separating the 
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two plasmas (resonant absorption). Ilonsonl (|1978l ) first utilized surface Alfven waves and resonant 
absorption as a mechanism to h eat coronal loo ps. T he transf e r of M HD wave energy by resonant 
absorption was also studied in iHollwea (jl9871 ) and IWentzell (I1979I). An alternativ e dissipation 
mechanism for surf ace Alfven waves is nonlinear wave steepening (jRudermanl Il992l ) . These and 
other efforts, e.g. (|Lee fc Roberts! Il986l ) have resulted in damping lengths which depend on the 
frequency of the waves, the nature of the magnetic interface, and the local plasma parameters 
(density, magnetic field and velocity). 

Utilizing these relations, the profile of the dampin g length in the wind has been estimated 
(jJatenco-Pereira Opherlll9 89; Narain &: Sharmalll9 98). All previous studies made assumptions 
about the wind. For example, (jNarain &: Sharmalll998l ) calculated nonlinear viscous damping of 
surface Alfven waves in polar coronal holes. They assumed two values of the super radial expansion 
of the magnetic field lines, profiles for density (based on observations), and a single frequency (0.01 
Hz). They obtained one profile, and concluded that the nonlinear damping of the surface Alfven 
waves in region of strong magnetic field expansion should contribute significantly to the heating in 
the solar wind. 

The surface Alfven wave damping length depends on the profile of the backgr ound Alfven 



speed . As was shown in a s urvey of Alfven speed profi les from several MHD models (lEvans et a 



2008 ). the Alfven profiles for Verdini Velli ( 2007 ) and Cranmer. van Ballegooijen fc Edgarl (|2007h 
are almost identical below 10 Rq, and the profiles were di f ferent from MHD models using empirical 
heating functions. The p rofile for lUsmanov &: Goldsteinl (|2006l ) is similar to these two, but differs 
very low in the corona. lEvans et al.1 (|2008l ) concluded that the inclusion of Alfven waves with 
empirical damping brought MHD models better in alignment with local heating studies that had 
the best agreement with observations. 

In the present study, we quantify the surface Alfven wave damping length for use in an MHD 
wave-driven model. We expect to find surface Alfven waves in the border between the fast and slow 
solar wind for two reasons: a) the gradient in density and b) the superradial expansion of the open 
magnetic field required to fill the space over closed streamers. We will focus on the superradial 
expansion of the field lines and compare our 3D MHD model with the observational study of 
Dobrzvcka et al.l (|1999l ). Our calculations show that waves with periods less than 1 hour (frequency 
greater than 2.8 xlO -4 Hz) are damped in the reg ion between 1-4 R q, the region of interest for 
both solar wind acceleration and coronal heating (jGrall et al.l ll996). No other major source of 
damping has been suggested for these waves in this region. We demonstrate the importance of 
the 3D geometry in our results. We show that the contribution of the damping of surface Alfven 
waves to the wind is on the order of magnitude as (or in some cases larger than) turbulence. It is 
important to note that this study of waves in a solar wind solution was not self consistent - we did 
not consider any back effects on the waves from the plasma. 



The paper is organized as follows: in section 2 we describe the theory and numerical simulation 
background. In section 3 we calculate the location of coronal hole boundaries, and the damping 



- 5 - 



length profile along those lines. In section 4 we estimate the energy surface Alfven waves will 
deposit below 4 Rq. We also calculate the heating invoked in a semi-empirical thermodynamical 
model, compare with the wave flux and discuss the results. Finally, conclusions can be found in 
section 5. 



Methods and Data 



2.1. Theory 



The inhomogeneity in density and/or magnetic field that gives rise to the surface Alfven wave s 
can be described analytically as either a discontinuity or finite layer (jHasegawa Uberoilll982l ). 
such as a flux tube. For the case of a rapidly expanding flux tube of width a, (where a is much 
smaller than the radius of the flux tube), surface Alfven waves form on the inner and outer surfaces. 
Thes e waves can interact . damp and deposit energy into the surrounding plasma with a damping 
rate (|Lee fc Robertslll986l ). 



sw 



ir(ka) 



8w 



(1) 



where k is the average wave number, u is the frequency and u\ and u>2 are the Alfven wave 
frequency on either side of the flux tube (1 representing inside and 2 outside.) We assume the 
width of the flux tube to be much smaller than the radius. This allows us to take ka = 0.1, as in 



Jatenco-Pereira fc Opherl (|1989l ) . If the frequency on the outside is much larger than the frequency 
inside (i.e., a strong inhomogeneity), then the damping rate is 



r 



sw 



7ruj(ka) 
4^/2 



(2) 



The surface Alfven wave damping length can be written as the Alfven speed va = y ^r; 



divided by the damping rate, 



L 



sw 



VA V A 



4V2 



Tsw 0J7r(ka) 
The initial damping length Lq can be written as 

VAO^Vl 



unr(ka) ' 



which, by taking ka = 0.1, can be simplified to 



L = 18 



VAO 



UJ 



(3) 



(4) 



(5) 



Utilizing the relation that a oc A(r)z tx r2 (where S is the superradial expansion factor of the 
field line), and fixing the frequency of the waves to be constant with height, the damping length in 
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the inertial frame is 



Lsw = Lq 



!*WJ* ] (1 + M A ) 



(6) 



where = ^^p^ is Alfven Mach Number, usw is the solar wind speed, B is the magnetic field 
strength, p is the mass density. The subscript indicates the variable is to be evaluated at the 
reference height. 

In the present study the model does not treat the chromosphere or the transition region; the 
lower corona is the inner boundary. The cell size at the inner boundary is jfg-R©, and so we chose 
ro=1.04i?Q as our reference height to assure that calculations would not include the solar surface 
inner boundary cells. 

We quantify the expansion of open field lines by S, given by. 

Sir) 



A cs (r) = A cs {r ) — 
Jo 



(7) 



where A cs (r) is the cross sectional area of the flux tube at distance r. A value of 2 for 5 indicates 
pure radial expansion. The lines which border closed field lines must open faster than radial to fill 
the space above the closed loops. In studi es where S is not a function of r, typical values chosen were 
2-6 for open field regions from 1-10 Rq (jMoore et al.lll99ll : iJatenco-Pereira. Opher fc Yamamoto 
1994 ; iNarain &: Sharmalll998l ). Here we determine S explicitly as a function of height from 1.04-10 

Rq. 

A similar parameter for the expansion of a field line is the superradial diverging factor or 
superradial enhancement factor /, as in 



A cs (r) = A cs (ro)(-) f(r). 



(8) 



We will use a 3D MHD model as a laboratory to estimate the contribution of the surface waves to 
the wind from 1.04-10 Rq. 



2.2. Generation of Steady State 



We obtain steady state solar wind solutions by using the Solar Corona compo nent of the 



Spac e Weather Modeling Framework (SWMF), developed by the University of Michigan (jToth et al 



2005). This 3D global magnetohydrodynamics (MHD) model incorporates Michelson Doppler Im- 
ager (MDI) magnetograms to g enerate an in i tial m agnetic field configuration with the Potential 
Field Source Surface model (see ICohen et al.1 (|2007l ) for details). An initial density is assumed on 
the solar surface (3.4 x 10 8 cm -3 ), and the MHD equations are evolved in local time steps (12,000 
iterations) and time accurate calculations (for ten minutes) to achieve steady state solutions for 
solar minima conditions in a 24 x 24 x 24 Rq domain. 
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The steady state was ge nerated with Carrington Rotation (CR) 1912. Solar wind solutions 
from SWMF were validated in lCohen et al.1 (|2008l ) from CR1916-1929 by comparing with Advanced 
Composition Explorer and Wind satellite data (near 1 AU). CR1912 w as chosen to allow for com- 
parisons to an observational study of the expansion of open field lines ([Dobrzvcka et al.lll999l ). 



Waves occur naturally as a perturbation to the MHD equations, and so their presence may be 
expected when solving the M HD equations in space and tim e. However, in global simulations waves 
have to be included explicitly (jUsmanov Goldsteinll2003l ) due to time and spatial limitations. The 
time step of this simulation (0.2 seconds) is less than the smallest period considered in this analysis 
(3 seconds). Additionally, the grid resolution is not enough to spatially resolve the waves. 



3. Coronal Hole Boundary Analysis 
3.1. Location and Expansion Factor 



Dobrzvcka et al.l (Il999l ) (herein referred to as D099) characterized the large-scale solar mag- 
netic topology during solar minima conditions (August 1996) using data from the Solar and Helio- 
spheric Observatory's (SOHO) Ultraviolet Coronagraph Spectrometer instrument. They analyzed 
the latitudinal dependence of two line emission intensities and found the values were constant 
within the large polar coronal holes but suddenly increased at the border of the holes and equa- 
torial streamers. D099 used this increase in intensity to identify the colatitude of coronal hole 
boundary (CHB), i.e. the border between open and closed magnetic field lines. In the present 
study, we identify the CHB locations as the open field line with the largest colatitude (the angle 
measured from the pole to the equator) in the steady state described in the previous section. We 
compare our calculation of the field expansion with the observed values from D099. 

Figure 1 shows the 3D coronal hole boundary (CHB) field lines obtained from the model. The 
different colors refer to: red as northeast (NE); blue as southeast (SE); green as southwest (SW); 
and purple as northwest (NW). The solar surface is shown colored by the radial component of the 
magnetic field. On 1996 August 17, there was an active region near the center of the disk, and 
mostly quiet sun at the intersection of the plane of the sky with the photosphere. 

The first two rows of Table [1] provide CHB colatitudes (angle measured from the pole to the 
CHB footpoint) found in D099 and this study. We find that 3 of the 4 simulated CHB are higher in 
latitude (i.e., smaller colatitude) compared to D099. As we discuss below, these differences could 
be due to projection effects. 

The superradial expansion factor S and superradial enhancement factor / were calculated as 
a function of height for each CHB field line. In Figure 2 we show / for the CHB lines in this study. 
Additionally, we have included the / profiles which correspond to the minimum and maximum 
asymptotic values of / for CHBs determined in D099. Rows 3 and 4 in Table Q] provide the 
asymptotic superradial enhancement factor value for each line. We find that both S and / cover 
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a larger range of values compared to DOB99. Only the SW line from our simulation falls in the 
range from D099 and only from 3-6 Rq. In general, we find that 3 of the 4 boundary lines (all 
except SW) have lower values in our simulation compared to D099. 

The best agreement for the location of the CHB between the studies is the SW line. The 
southern hemisphere from the simulation is more similar to the results deduced f rom observations 



i n the plane of the sky in D099 than the northern hemisphere. Overall, we find that lDobrzvcka et al 



(|1999l ) had larger values of /, and a small range of values, for the same field lines. As we show 



below, this result is due to 3D vs. 2D projection effects. 



3.2. Projection Effects 

In order to quantify the significance of 2D projection effects in the context of the expansion of 
field lines, we calculated / for four field lines shown in Figure 3a. The projection of the field lines 
on the plane of the sky is shown in Figure 3b. Contours of solar wind speed (in ^p) are shown in 
Figures 3a and b, and the solar surface is shown as white sphere. Figure 3c provides the superradial 
enhancement factor /, calculated according to Equation [8] for each field line (labeled A, B, C and 
D as in Figure 3b). The calculation of the 3D line is shown as solid lines, and the 2D projection 
line calculation is shown as dashed lines. 

Figure 3c demonstrates that a 2D projection of a 3D field line on the plane of the sky can 
overestimate the divergence of the line. The 3D lines from A-D all approach values of f 4 while 
the 2D estimates vary between 4 and 13. This projection effect explains why the values for / from 
our simulation are smaller than those determined in the observational study of D099. It is crucial, 
therefore, to do studies of surface Alfven wave damping in a 3D simulation in order to capture the 
true divergence of the field lines. 



3.3. Damping Length 

In Figure 4, we plot L$w (see Eq. [6]) which was calculated using parameters p, B and u sw 
from the steady state solution. Figure 4a presents L$w f° r the coronal hole boundary field lines 
in Figure 1 with frequency 4.17 x 10 -3 Hz, normalized to the initial damping length Lq of the SW 
line (chosen because of the agreement with D099). This normalization allows for comparison of 
the profile features from different source regions as a function of height. In Figure 4b we feature 
only the SW line and present L$w f° r several frequencies, from 3.3 x 1CP 1 — 3.8 x 1CP 6 Hz. It 
can be seen in Figure 4b that frequencies above 2.8 x 10 -4 (short dashed line) will be appreciably 
damped within a few solar radii of the surface. 

Figure 4a shows distinctly different profiles from the southern and northern CHB lines. We 
examined the source region of each footpoint and found that the SE and SW lines originated near 



small active regions in which the radial component of the magnetic field B r ~ 50 G. Both northern 
hemisphere li nes o riginated from quiet sun regions (B r ~ 1 G). For SWMF and other MHD models, 



Evans et al.l (|2008l ) showed that the Alfven profile will contain a maximum, or hump, if the source 
region is quiet sun. The profile from an active region in global models begins at a maximum value, 
and drops to less than a few hundred within one solar radius from the surface. 

The profile of Lsw is controlled by the Alfven speed profile. The normalized profiles in Fig. 
4a show that the position corresponding to Lsw = 1-R© is closest to the Sun for active regions. 
The profiles from quiet Sun source regions have a plateau, pushing Lsw = l-R© further from the 
Sun. The implication of this result can be seen in the equation relating the Alfven wave energy 
density, 



+ Mad \ 2 } 

+ M A ) 6XP V L sw 



If the damping length is IRq or less, then the waves will be damped close to the Sun. Therefore, 
the presence of the hump means the energy of the surface Alfven wave can travel further into the 
corona before substantial damping occurs. This means that the quiet sun region will damp more 
surface waves at further distances, so it is more efficient in carrying the wave momentum out into 
the corona. Active regions will damp closest to the Sun. 



4. Dissipation of Wave Energy and Heating 
4.1. Wave Energy 



In the previous section we considered h ow surface Alfven waves at the sol ar surface with the 
frequency range 3.8 x 10 -6 to 3.3 x 10 _1 Hz (jCranmer h van Ballegooijenll2005l ) would be damped 
in our background solar wind environment. We found that waves with frequency above 2.8 x 10~ 4 
Hz were appreciably damped below 4i? Q . We now consider the contribution of their wave flux to 
the energy of the wind. There will be a contribution to the momentum of the wind as well, but in 
this analysis we ignore this contribution. We assume that the wave damping will contribute solely 
to heating the wind; therefore we derive here an upper limit on their contribution to the heating o f 
the plasma. The spectra of the surface Alfven waves (IJatenco-Pereira. Opher &: Yamamotdll994l ) 
is: 

<t> AW {uj) = 00 (~) ~ a T*L (10) 



where 



1.3 x 10 5 ^ 



\ujJ cm 2 sHz 

rm ^sHz ' ^ 1S ^ ne mean frequency in the observed range and the power index 



corresponding to the low frequency waves we are considering is a=0.6 (|Tu et al.lll989l ) 



We assume that this flux of Surface Alfven waves is propagating along open field lines during 
solar minima. This flux will decrease with distance as 
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where Lgw is the damping length, for each frequency. With the damping lengths calculated in 
the previous section, we calculate how much flux is lost between R = 1.04 (our self-imposed 
reference height) and 4i? . The height of 4i? was chosen because the contribution to solar wind 



acceleration and coron al heating must be deposited within a few solar radii of the Sun (|Kohl et al 



19981 : lEsser et aUll99a ). 

The contribution to the solar wind at each frequency is 



<t>lost(u) = 00 



1 — exp 



V 



12 



-dr 



<>(«)(*)*( 



V A 
VAO 



;i + M A ) 



(12) 



where the limits are n = 1.04i?Q and r<i = 4i? Q , and the definition of L sw from Equation [6] has 
been included (va, S, Ma, are all functions of r.) 



We replace Lq i n Eq. with Eq. [SJ The total flux lost is found using 



'lost.total 



\L0 



1 — exp 




-UJ- 



s 

nvAo 



-dr 



18r*v 2 A (l+M A ) 



duj 



(13) 



where the limits are u\ = 2.8 x 10~ 4 Hz and uo 2 = 0.3 Hz. 

Next we compare the potential contribution of the wave flux to the heating in the model (Eq. 



22|) . It should be stressed that we are not doing a self-consistent calculation: we are estimating 
wave flux from a model that does not include waves, and we are not considering any feedback of 
the waves on the plasma. 



4.2. Heating of the Corona 



The simulation analyzed in this paper is that of a thermally driven solar wind. The polytropic 
index V in the model varies in space by utilizing the Wang-Sheeley-Arge model and the conservation 
of energy along a solar wind field line using the Bernoulli equatio n. This serves to arti ficially heat 
the wind ([ Cohen et al.l 120071 ) in a manner mimicking turbulence iRoussev et al.l (|2003l ) . Figure 5 
shows the distribution of V in the plane of the sky on 1996 August 17. In this section we quantify 
the heating due to this variable gamma, and do a comparison with the energy deposited by damped 
surface Alfven waves, as in Eq. [131 

The first law of thermodynamics attributes changes in the internal energy U of a gas to work 
done on or by the gas W, and heat added to or removed from the gas Q. In the case of an ideal 
gas, the change in internal energy can be written as dU = c v dT, where c v is the specific heat at 
constant volume. The work is expressed as dW = —pdV, where p is the pressure and V is the 
volume. The first law can therefore be written as 



dQ = c v dT + pdV. 



(14) 
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By introducing the ideal gas equation of state and assuming that the ratios of specific heats are 
constants, one can derive a polytropic equation, 

p p 

C P~ C _ pa 
pc v —c r 



const. 



(15) 



where a is referred to as the polytropic index. The notation stems from iParkerl (|1963l ) to clarify 
that this index can (but need not) be the ratio of specific heats, and that we are not necessarily 
considering an adiabatic process. The symbol 7 is typically used for the ratio of specific heats, 
and in the case of an adiabatic expansion (no heating enters or leaves the system), a = |. An 
isothermal wind expansion would be characteri zed by a = 1. Observations of the solar wind have 



indicated that a=1.46-1.58 in the heliosphere (jTotten. Freeman fc Arvalll995l ). A value closer to 
unity is adopted in some global MHD models in the region near th e Sun in order to generate fa st 
solar wind and match temperature observations in the heliosphere (jUsmanov &: Goldstein! 120031 ). 



All previous discussion had the underlying assumption that a was constant with height. If that 
condition is not met, the n the polytropic index is referr ed to as an effective (or local) polytropic 
index and written as T (jTotten. Freeman &: Arval Il995l ). The polytropic equation (Eq. [To"]) is 
modified to 

^dlnp 



dlnP 



+ Inp 



dT 



(16) 



dr dr dr 

such that the relationship between density and pressure is not simple. The variation of V with height 
has been utilized to heat the solar wind used in this paper. We will characterize the additional 
heating provided by the prescribed distribution of V in our solar wind simulation, and argue that 
surface Alfven waves damped near the Sun could replace this artificial heating and move the model 
towards a more physical treatment of the solar environment. For a solar wind with a cons tant ratio 
of specific heats 7, the conservation of energy can be written as (jManchester et al.1 12004 ): 



de 
dt 



+ V 



u 



B 2 \ 



u 



B)B 



47T 



(17) 



where p is the thermal pressure, q is the additional heating function, and the energy density is 

pu 2 v B 2 



+ 



7 — 1 57T 

Recent global MHD studies adopt an exponential function for the form of q with several free 
parameters in ord e r to benchmark the mo del with observations during solar minima conditions 
(|Groth et al-lboool ; iManchester et all I2OO4I ) . Substituting Eq. [18] into Eq. [TTl and setting time 
derivatives to zero for a steady solar wind, we find: 

pu 2 B 2 \ fu B)B 



V- 



u 



IP 



B 2 \ 



pg-u + - 



(19) 



^ 7 - 1 2 4"7T / 47T 

As discussed at the beginning of this section, there is no heating function q in the simulation used 
in this paper, and the ratio of specific heats 7 is replaced by the effective gamma T, 



V ■ 



u 



Tp 



+ 



pu 



B 2 \ 
+ 4Tj 



u 



B)B 



47T 



pg-u. 



(20) 
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Although r has both latitudinal and azimuthal dependence below 4Mq, we consider only at the 
radial variation, and so we replace V by We assume that we have exactly the same solar 
wind solution in the two cases we are considering: that of a variable gamma and of an additional 
volumetric heating function with 7 = |. In order to quantify the amount of heating in the model 
with variable gamma, we subtract equation [20] from [T9j 



d_ 

dr 



Tp 



IP 



r- 1 7-1 



(21) 



Equation [21] can be written as: 



d (u r p) 
dr 



r-i 7-1 



dr u r p 

dr (r - 1) 



(22) 



Knowing how T, p and u r vary along any radial line, and setting 7 = | we can integrate 
equation [22] between 77=1. 04R and r2=4R to find the heat input along any field line: 

Q=[ r2 ldr^- s . (23) 

This equation gives the heat deposited into the system between the two heights. We compare Q 
with the flux of damped surface Alfven waves (Equation 1131) 

Table [T] provides (friost an d Q for the CHB field lines. The expansion / of the NW line from the 
simulation has the best match to the observations (5.55 compared to 6.0). The surface Alfven wave 
flux along this line, and along the NE lin, is larger than the heating Q by an order of magnitude. 
The geometrical properties of the SW line also match well with observations, however the wave flux 
for it (and also the SE line) account for 20 % of the required heating. The distinction between the 
southern and northern lines is the source region: they come from a stronger magnetic field region. 
Near an active region, the second term in Eq. [22] (which includes the pressure and radial velocity) 
is larger than a quiet sun region. Therefore, we expect Q to be larger than the surface Alfven wave 
flux along lines from active regions. 

As we assumed all of the wave flux goes to heating, this procedure gives an upper limit on 
the contribution of the damping of surface Alfven waves along an open magnetic field line to the 
heating along that line. A random sampling of 7 open field lines in the northern hemisphere with 
footpoints in the plane of the sky (see Figure 3) yielded <pi os t that were on the order, or an order 
of magnitude larger than the Q. 



5. Conclusions 



This work is the first study to look at surface Alfven waves in a 3D environment without 
assuming a priori a geometry of the field lines or magnetic and density profiles and strengths. We 
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showed the calculation of the expansion of field lines must be done in a 3D environment. Our 
calculations show that waves with periods less than 1 hour (frequency greater than 2.8 xlO -4 Hz) 
are damped in the re gion between 1-4 Rq, the region of interest for both solar wind acceleration 



and coronal heating (jGrall et al.lll996l ). We showed that the quiet sun region will damp surface 
waves at further distances, so it is more efficient in carrying the wave momentum out into the 
corona. Surface waves formed on flux tubes with footpoints in an active region will damp closer to 
the Sun. The required heating from an active region was found to be larger than the damping of 
surface Alfven wave flux, therefore another mechanism (such as turbulence) may be the dominant 
heating in these regions, with surface Alfven waves contributing approximately 20% of the heating. 

We estimated damping of surface Alfven waves in the border between open and closed field 
lines at heights 1.04-4 Rq due to the superradial expansion of the field lines. As some of the wave 
flux would go to the momentum of the wind, we provide an upper limit on the contribution of 
surface Alfven waves to the heating of the solar wind. In the region between open and closed field 
lines, within a few solar radii of the surface, no other major source of damping has been suggested 
for the low frequency waves we consider here. 

Our results demonstrate that it is not necessary to have turbulence in order to heat the solar 
wind - and that it is imperative to include the physics of surface Alfven wave damping in solar 
wind models in order to more physically model the heating. Surface Alfven waves cou ld also be 



pr 
2C 



rese nt in the solar wind, in the flux tube structures said to fill interplanetary space (jBorovsky 



2008). Another environment which could support these waves are Coro tating Interactio n Regions 



(CIRs), due to the inhomogeneity in density present in these structures (|Tsubouchill2009l ) . Both of 
these topics will be addressed in future works. 

It is important to note that this study of waves in a solar wind solution was not self consistent 
- we did not consider any back effects on the waves from the plasma. We simply tried to estimate 
if the waves could produce the heating required to create the solar wind solution from the model. 
In the future we will pursue other damping mechanisms with the goal of incorporating the key 
mechanisms of wave damping in self-consistently in global MHD models to improve the lower 
corona. 



We would like to thank NASA Ames for the use of the Columbia supercomputer. This research 
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Fig. 1. — 3D coronal hole boundary field lines obtained from the simulation for 1996 August 17. 
The different colors refer to: red as northeast line (NE); blue as southeast line (SE); green as 
southwest line (SW);and purple as northwest line (NW). The black lines show closed equatorial 
streamers, and the arrows give the direction of the magnetic field. The solar surface is shown 
colored by the radial component of the magnetic field. 
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Fig. 2. — The superradial expansion factor / (see Eq. [8]) for the lines in Fig. [T] (colors corre- 
spond) from 1.04-10 solar radii (Rs). We also show the profiles which correspond to the minimum 
(dashed line) and max imum (dash-dot line) values of / determined in the observational study by 
Dobrzycka et al.l (119991 ). We find that / covers a larger range of values compared to the observa- 
tional study. 
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Fig. 3. — Four field lines (A, B, C and D) shown (a) in 3D and (b) as their projections on the 
2D plane of the sky. Contours of solar wind speed (in ^p) are shown on the plane, and the solar 
surface is shown as a white sphere. The expansion factor calculated from the 3D field line and 2D 
projections for each line is shown in (c). The 2D projection can overestimate the expansion factor 
as compared to the 3D calculation. 
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Fig. 4. — a) Surface Alfven damping length (L sw ) profiles for the coronal hole boundaries in Fig. 
1 (colors correspond), normalized to the Lq of the SW line (green) from 1.04-10 solar radii (Rs). 
Note the plateau in the profiles for lines whose source region on the sun is quiet sun, differing from 
field lines with footpoints near active regions, whose profiles drop quickly, b) Profiles for the SW 
coronal hole boundary field line corresponding to different frequencies. From the bottom up to top 
profile: 3.3 x 10 _1 Hz (dash-dot-dot line); 1 x 10~ 2 Hz (dash-dot line); 4.17 x 10 -3 Hz (dot line); 
1.67 x 10~ 4 Hz (long dash line); 2.8 x 10~ 4 Hz (short dash line); 1.5 x 10~ 5 Hz (thin solid line); 
3.8 x 10~ 6 Hz (thick solid line). Waves with frequencies above 2.8 x 10 -4 are appreciably damped 
below 4Rs. 
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Fig. 5. — The distribution of T, the effective polytropic index, in the plane of the sky on 1996 
August 17. r is specified on the solar surface (shown as white sphere) using the Bernoulli integral, 
and is interpolated to 1.1 at 2.5 Rs (inner black circle). Above 2.5Rs, T varies linearly until 12.5Rs 
(outer black circle), above which it has the value 1.5. The white lines show the boundaries of grid 
refinement. 



-20- 



Jatenco-Pereira, V. k Opher, R. 1989, A&A, 209, 327 
Jatenco-Pereira, V., Opher, R. k Yamamoto, L. C. 1994, ApJ, 432, 409 

Jess, D. B, Mathioudakis, M., Erdelyi, R., Crockett, P. J., Keenan, F. P. k Christian, D. J. 2009, 
Science, 323, 1582 

Kohl, J. L. et al. 1998, ApJ, 501, L127 

Leake, J. E., Arber, T. D. k Khodachenko, M. L. 2005, A&A, 442, 1091 
Lee, M. A. k Roberts, B. 1986, ApJ, 301, 430 
Lionello, R., Linker, J. A. k Mikic, Z. 2009, ApJ, 690, 902 
Manchester, W. B. et al. 2004, J. Geophys. Res., 109, A02107 

Matthaeus, W. H., Zank, G. P., Oughton, S., Mullan, D. J. k Dmitruk, P. 1999, ApJ, 523, L94 

Matthaeus, W. H., Dmitruk, P., Oughton, S. k Mullan, D. 2003, SOLAR WIND TEN: Proceedings 
of the Tenth International Solar Wind Conference, 10.1063/1.1618627 

Mikic, Z., Linker, J. A., Schnack, D. D., Lionello, R. k Tarditi, A. 1999, Physics of Plasmas, 6, 5 

Moore, R. L., Musielak, Z. E., Suees, S. T & An, C. -H. 1991, ApJ, 378, 347 

Narain, U. k Sharma, R. K. 1998, Solar Physics, 181, 287 

Ofman, L. k Davila, J. M. 1997, J. Geophys. Res., 103, A10, 23677 

Okamoto, T. J. et al. 2007, Science, 318, 1577 

Parker, E. N. 1963, Interplanetary Dynamical Processes, Wiley-Inter science, New York 

Parker, E. N. 1965, Space Sci. Rev., 4, 666 

Parker, E. N. 1991, ApJ, 376, 355 

Patsourakos, S. k Klimchuk J. A. 2009, ApJ, 696, 760 

Roussev, I. I. et al. 2003, ApJ, 595, L57 

Ruderman,M. S. 1992, ApJ, 399, 724 

Suzuki, T. K. k Inutsuka, S. 2005, ApJ, 632, L49 

Suzuki, T. K. k Inutsuka, S. 2006, J. Geophys. Res., Ill, A06101 

Tomczyk, S. et al. 2007, Science, 317, 1192 



- 21 - 



Tomczyk, S. & Mcintosh, S. W. 2009, ApJ, in press 
Toth, G. et al. 2005, J. Geophys. Res., 110, A12226 

Totten, T. L., Freeman, J. W. & Arya, S. 1995, J. Geophys. Res., 100, Al, 13 
Tsubouchi, K. 2009, J. Geophys. Res., 114, A02101 

Tu, C. Y., Marsch, E. & Thieme, K. M. 1989, Geophys. Res. Lett. 17, 283 
Usmanov, A. V. & Goldstein, M. L. 2003, J. Geophys. Res., 108, A9, 1354 
Usmanov, A. V. & Goldstein, M. L. 2006, J. Geophys. Res., Ill, A07101 
Verdini, A. & Velli, M. 2007, ApJ, 662, 669 
Wentzel, D. G. 1979, ApJ, 227, 319 
Wentzel, D. G. 1989, ApJ, 336, 1073 



This preprint was prepared with the AAS IATjrjX macros v5.2. 



-22- 



Table 1: Properties of Coronal Hole Boundary Field Lines 





NE 


NW 


SE 


SW 


6>n a (Dobrzvcka et al. 1999) 


29.0° 


31.0° 


23.3° 


28.0° 


#o This Study 


37.2° 


22.3° 


19.2° 


25.2° 


f,n B . b (Dobrzvcka et al. 1999) 


6.56 


6.00 


7.30 


6.5 


f WR(3 This Study 


2.56 


5.55 


4.71 


7.62 


Q c (^-) 

^ v cm z s > 


8.9 x 10 3 


8.0 x 10 3 


2.7 x 10 5 


3.1 x 10 5 


4>lost ( cm 2 s ) 


6.1 x 10 4 


6.2 x 10 4 


6.4 x 10 4 


6.4 x 10 4 



"Colatitude of the line - angle measured from the pole to the footpoint on the solar surface. 
'Value for the superradial enhancement factor at R=10 Rq 
c Heating calculated along the field line (see section 4.2) 
d Alfven wave flux deposited into the wind (see section 4.1) 



